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ABSTRACT 

Using 3.7 years of Fermi-L AT data, we examine the diffuse 80 — 200 GeV emission in the inner 
Galaxy and find a resolved gamma-ray feature at ~ 110 — 140 GeV. We model the spatial distribution 
of this emission with a ~ 3° FWHM Gaussian, finding a best fit position 1.5° West of the Galactic 
Center. Even better fits are obtained for off-center Einasto and power-law profiles, which are preferred 
over the null (no line) hypothesis by 6.5cr (5.0a/5.4cr after trials factor correction for one/two line case) 
assuming an NFW density profile centered at (£,b) = (— 1.5°, 0°) with a power index a = 1.2 . The 
energy spectrum of this structure is consistent with a single spectral line (at energy 127.0 ± 2.0 GeV 
with x 2 = 4.48 for 4 d.o.f.). A pair of lines at 110.8±4.4 GeV and 128.8±2.7 GeV provides a marginally 
better fit (with x 2 = 1-25 for 2 d.o.f.). The total luminosity of the structure is (3.2 ±0.6) x 10 35 erg/s, 
or (1.7 ± 0.4) x 10 36 photons/sec. The energies in the two-line case are compatible with a 127.3 ± 2.7 
GeV WIMP annihilating through 77 and 7Z (with x 2 — 1-67 for 3 d.o.f.). We describe a possible 
change to the Fermi scan strategy that would accumulate S/N on spectral lines in the Galactic center 
4 times as fast as the current survey strategy. 

Subject headings: gamma rays — diffuse emission — milky way — dark matter 



1. INTRODUCTION 



Although various cosmological and astrophysical ob- 
servations provide compelling evidence for dark matter 
(DM), which constitutes ~ 80% of the matter in the Uni- 
verse, we still know little about its intriguing n ature (e.g. 



iBertone et al.ll2005 



ittie about its intriguing n ature (e.g. 
Hooper fe Profumoll2007 ). Among a 



forest of dark matter models, stable Weakly Interacting 
Massive Particles (WIMP) have been predicted in many 
extensions of the Standard Model of particle physics (e.g. 
BergstromlboOOh . The WIMP with the virtues of weak 
scale masses and couplings is an excellent dark matter 
particle candidate whi ch can annihilate into high en- 
ergy gamma-rays (e.g. IBergstrom et al.lll998[ ). The in- 
ner Galaxy provides one of the most promising regions 
on the sky t o search for WIMP annihilation produced 



gamma rays (Ackerm ann et al.l 12011c lAbramowski et al 



20111 : iGondolo fc Silkl Il999h . The expected relatively 
higher annihilation rate due to higher dark matter par- 
ticle density provides a potential window to identify any 
non-gravitational dark matter signatures. 

The "smoking-gun" signal of annihilating dark mat- 
ter would be a monochromatic gamma-ray line (or lines) 
in a region of high dark matter density, either in lo- 
cal dwarf galaxies or in the Galactic center. This line 
could be produced by dark matter decays or annihila- 
tions into two photons, or two-body final states involv- 
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ing one photon plus a Higgs boson, Z boson, or other 
chargeless non-SM particle. In most models, dark mat- 
ter does not annihilate directly to photons, but in models 
where it annihilates to charged lepton pairs, there may 
be loop interactions that produce photons in two-body 
final states, yielding spectral lines. These lines can pro- 
vide a signature of dark matter even in a complicated 
astrophysical environment, because no known astrophys- 
ical process produces gamma-ray lines at E ^> 1 GeV. 
However, these loop processes would be suppressed by 
1-4 orders of magnitude compared to the total a nnihi- 
lation or decay rate (e.g. IBergstrom fc Ulliolll997l ). The 



unprecedented sensitivity an d spectral resolution of the 
Large Area Telesco pe (LAT: iGehrels fc Michelsonlll999 : 



Atwood et al.ll2009f ) aboard the Fermi Gamma-ray Space 
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Telescope make it possible to search for dark matter an- 
nihilation lines over the whole sky up to a few hundred 
GeV. 

However, other signals can mimic a weak line. Features 
such as a spectral edge or a broken power law can be mis- 
taken for a line in noisy data, especially when smoothed 
by the instrumental response. Thus careful separation of 
possibly diffuse dark matter emission from other diffuse 
components is crucial. 

For example, the Fermi bubbles exten d ~10 kpc abov e 
and below the Galactic center (GC) ( Su et all |2010| ). 
Their emission has a flat energy spectrum in E 2 dN/dE 
from hundreds of MeV to ~100 GeV implying the abil- 
ity of the Fermi bubbles to accelerate cosmic ray (CR) 
electrons up to ~TeV, if the gamma rays are produced 
by inverse Compton (IC) scattering of the Cosmic M i- 
crowave Background (CMB) photons (|Su et alJlioiol ). 
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Fig. 1. — All-sky Fermi-LAT 3.7 year sky maps in 4 energy bins ranging from 100 to 180 GeV. We use CLEAN event class and point sources 
have been subtracted based on the Second Fermi-LAT catalog (2FGL). Large sources, including the inner disk (—2° < b < 2°,— 180° < 
£ < 180°), have been masked. The maps have been smoothed for display with a Gaussian kernel of FWHM = 2°. 
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Fig. 2. — Same as Figure [TJ but with SOURCE class events. This event class contains substantially more background. 



If there is a spectral break in the Fermi bubble spec- 
trum at high energy, the superposition of such a spec- 
trum with other soft-spectrum diffuse components might 
mim ic a bump feature on t op of a continuum power 
law (jProfumo fc Linden! 120121 ). 



In this work, we use 3.7 years data of LAT to study 
the diffuse gamma-ray emission toward the inner Galaxy 
at 80 < E 1 < 200 GeV. In Section [2] we describe our 
LAT data selection and map making. In Section O we 
show that the gamma-ray maps reveal a novel gamma- 
ray cusp toward the Galactic center. We characterize the 



morphology of the gamma-ray cusp and employ regres- 
sion template fitting to determine its energy spectrum 
in Section [U In Section [5] we study the detailed spatial 
profile of the cusp structure by fitting various templates. 
Section [6] contains analysis of the spectral line profile, 
and searches for instrumental effects in photon samples 
from the Galactic plane and the Earth's limb. We de- 
rive the energy spectrum of the cusp assuming various 
dark matter density profiles in Section [3 We propose a 
modified survey strategy in Section [8] that would allow 
Fermi-LAT to confirm this gamma-ray line signal in ~1 
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Fig. 3. — All-sky CLEAN 3.7 year maps in 5 energy bins, and a residual map (lower right). The residual map is the 120 — 140 GeV map 
minus a background estimate, taken to be the average of the other 4 maps where the average is computed in E 2 dN/dE units. This simple 
background estimate is sufficient to remove the Galactic plane and most of the large-scale diffuse structures and even bright point sources. 
A cuspy structure toward the Galactic center is revealed as the only significant structure in the residual gamma-ray map. All of the maps 
are smoothed with a Gaussian kernel of FWHM = 10° without source subtraction. 

ment response functions, known as Pass 7 (P7_V6) 4 . For 



year with no trials factor, and we summarize our main 
findings in Section [9l 

2. MAP CONSTRUCTION 
For this project, we constructed full-sky maps from the 



LAT event files as in our previous work ([Pooler et al 



2010l : ISu et aHl2Qld : ISu fc Finkbeinerll2Q12h . except that 



we now use 3.7 years of Pass 7 (P7_V6) data. 

2.1. Fermi data selection 

The Fermi LAT is a pair-conversion telescope, in which 
incoming photons convert to e + e _ pairs, which are then 
tracked through the detector. The arrival direction and 
energy of each event are reconstructed, and the time of 
arrival recorded. Event files for every week of the mission 
are available on the Internet, and it is from these files 
that we build our maps. 

The point spread function (PSF) is about 0.8° for 68% 
containment at 1 GeV and decreases with energy as r^g ~ 
E~ - 8 , asymptoting to ~ 0.2° at high energy. The LAT 
is designed to survey the gamma-ray sky in the energy 
range from about 20 MeV to several hundreds of GeV. 

We use the latest publicly available data and instru- 



most figures in this work we use the CLEAN event class, 
which has larger effective area than ULTRACLEAN and 
lower background than SOURCE. In a few cases, we show 
figures made with ULTRACLEAN or SOURCE events as ev- 
idence that this choice has no qualitative effect on our 
results. 

Photons coming from the bright limb at Earth's 
horizon, dominantly produced by grazing-incidence CR 
showers in the atmosphere, are a potential source of con- 
tamination. We minimize this background by selecting 
events with zenith angle less than 100° as suggested in 
the Fermi Cicerone 5 . We also exclude some time in- 
tervals, primarily while Fermi passes through the South 
Atlantic Anomaly. 

2.2. Map making 

We generate full-sky maps of counts and exposure us- 
ing HEALPix, a convenient equal-area iso-latitude full- 

4 Details at http://fermi.gsfc.nasa.gov/ssc/data/analysis/ 
documentation/Pass7_usage . html 

5 http: //fermi .gsfc .nasa.gov/ssc/data/analysis/documentation/. 
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Fermi 120 < E < 140 GeV 



Fermi 140< E< 160 GeV 




Fig. 4. — The same as Figure [3] but we have subtracted point sources before smooth the maps. 



sky pixelization widely used in the CMB community. 6 
Spherical harmonic smoothing is straightforward in this 
pixelization, and we smooth each map by the kernel re- 
quired to obtain an approximately Gaussian PSF of some 
target FWHM, usually 10°. We generate maps for front- 
and back-converting events separately, smooth them to 
a common PSF, and then combine them. 

We construct maps both with and without point source 
subtraction. We subtract point sources listed in the Sec- 
ond Fermi-L AT catalog (2FGL), which is based on 24 
months of P7_V6 LAT observations. 7 The PSF and ef- 
fective area of the Fermi-LAT varies with energy, and we 
subtract each point source from the maps in each energy 
bin, using the in-flight version of the PSF contained in 
the P7_V6 IRFs. 

For the 400 brightest and 400 most variable sources, 
the subtraction is noticeably imperfect at lower energies 
(and we assume it is also at the higher energies used in 
this work), so we interpolate over the core of the PSF 

6 HEALPix software and documentation can be found at 
http://healpix.jpl.nasa.gov, and the IDL routines used in 
this analysis are available as part of the IDLUTILS product at 
http : / / sdss3data . lbl . gov/ software/ idlutils. 

7 http: //fermi .gsfc .nasa.gov/ ssc/data/access/lat/2yr_catalog, 
the file we used is gll_psc_v07 . f it 



after subtracting the best estimate. We also mask out 
sources including Geminga, 3C 454.3, and LAT PSR 
J1836+5925 and large sources like Orion and t he Magel- 



lanic Clouds as we did in th e previous papers ([Su et al 



201(1 ISu fc Finkbeinerll2Q12f ) although they are unlikely 
to be a problem at E > 80 GeV, where we are searching 
for lines. 

We produce the exposure maps using the gt It cube 
and gtexpcube2 tasks in the Fermi Science Tools. For 
bright sources, the exposure is set to zero for excised 
pixels. For the smoothed maps, both the count map 
and exposure map are smoothed, and then divided. At 
high energies, where the PSF is small, this effectively 
interpolates over the masked pixels. 

3. A MODEL-INDEPENDENT SEARCH FOR 
LINE-EMITTING REGIONS 

We begin by constructing maps in broad (20 GeV) 
energy bins and taking linear combinations that cancel 
most Galactic emission. This has the potential to re- 
veal a component with an unusual spectrum in a model- 
independent way. 

In Figures Q] and [2j we show the full sky map in four 
energy bins in the range 100 to 180 GeV, using CLEAN 
and SOURCE events, respectively. Even after nearly four 
years of observation, the number of gamma-ray photons 
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Fig. 5. — Residual maps the same as the lower right panel of Figures [3] and [4] but using SOURCE events before (upper left panel) and 
after (upper right panel) point source subtraction, or using ULTRACLEAN events before (lower left panel) and after (lower right panel) point 
source subtraction. 
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Fig. 6. — Residual maps by subtracting the average map of 
100 — 120 GeV and 140 — 160 GeV maps from the average map of 
80 — 100 GeV and 160 — 180 GeV maps (individual maps are shown 
in Figure [3} • Most of the large scale diffuse gamma-ray structures 
(including the Galactic plane) and bright point sources have been 
removed, and no gamma-ray cuspy structure toward the Galactic 
center is visible in the residual gamma-ray maps. The maps are 
constructed using Fermi-L AT 3.7 year CLEAN events (upper panel) 
and SOURCE events (lower panel), respectively. 

with E > 100 GeV is still quite limited and the maps 
are Poisson noise dominated. In order to inspect diffuse 
gamma-ray structure, we smooth the maps with a Gaus- 
sian kernel of FWHM = 10°. Smoothed CLEAN maps in 



Fig. 7. — Residual map after subtracting the 120 - 140 GeV 
ULTRACLEAN map from the SOURCE map, with both smoothed to 
FWHM = 10° . The map is consistent with Poisson noise, with no 
excess toward the Galactic center. 

various energy bins with E > 80 GeV, before and after 
subtracting point sources, are shown in Figures [3] and [4j 
respectively. To look for any diffuse gamma-ray emis- 
sion component with a distinctive energy spectrum, we 
examine various linear combinations of maps that cancel 
out the Galactic plane emission and visually inspect the 
residual maps. 

Interestingly, when we subtract the average map (aver- 
aged in E 2 dN/dE units) of the 80 - 100 GeV, 100 - 120 
GeV, 140 - 160 GeV, and 160 - 180 GeV maps from 
the 120 — 140 GeV map, most of the large scale diffuse 
gamma-ray structures (including the relatively bright 
Galactic plane) and visible point sources have been 
largely removed. However, a resolved cuspy structure to- 
ward the inner Galaxy with a size slightly larger than the 
smoothing kernel is revealed as the only visible structure 
in the residual gamma-ray map. We have repeated the 
analysis only subtracting 80 - 100 GeV and 100 - 120 
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Fig. 8. — Spatial templates used in the Poisson likelihood analysis. Upper left: Galactic disk template, upper right: Fermi bubble 
template, lower left: gamma-ray cuspy template as a Gaussian distribution with FWHM = 4°, lower right: outer ring template as a 
Gaussian distribution with FWHM = 10°, but masking out the central region where the gamma-ray cuspy template is. In Figure [9] and 
Figure [TTI we split the Fermi bubble template into two components one with \b\ > 30° and the other with \b\ < 30°. 



GeV maps from the 120 — 140 GeV map, or only sub- 
tracting 140 - 160 GeV and 160 - 180 GeV maps, and 
we obtained similar residual structure toward the Galac- 
tic center. We also repeat this exercise using SOURCE 
and ULTRACLEAN event classes and the residual maps are 
shown in Figure [5l The GC excess is similar in each 
case, and the rest of the gamma-ray sky is consistent 
with Poisson noise. This excess at ~ 120 — 140 GeV 
strongly suggests a novel diffuse gamma-ray component 
toward the Galactic center with unusual spectrum. Fur- 
thermore, this energy range coincides with the recently 
sugge sted tentative signature of gamma-ray excess at 130 
GeV (Bringm ann et al 1l2012l : IWenigerll2012h. which is un- 



der active debate in the literature (|Bovarsky et al.ll2012 



Tempel et allboii IProfumo fc Linden! 12 0121 ). Assuming 
the distance to the Galactic center is Rq = 8.5 kpc, the 
size of the gamma-ray cusp is < 1 kpc. 
We show in Figure [6] the difference map between the 



average map of 80 — 100 GeV and 160 — 180 GeV maps 
and the average map of 100 - 120 GeV and 140 - 160 
GeV maps. The difference map is consistent with Pois- 
son noise and no diffuse gamma-ray excess toward the 
inner Galaxy is visible. In order to test whether the ex- 
cess is due to residual cosmic ray contamination, we sub- 
tract ULTRACLEAN sky maps from SOURCE sky maps. This 
residual map should be mostly dominated by cosmic rays 
since a large fraction of the real gamma-ray photons have 
been removed. Indeed, Figure [3 demonstrates that there 
is no excess toward the inner Galaxy in this map, thus 
we can rule out the possibility that the central excess is 
due to cosmic ray contamination in the LAT data. 

Toward the inner Galaxy, the Fermi bubbles extend 
~ 50° above and below the Galactic center, with a width 
of ~ 40° in longitude. The gamma-ray emission associ- 
ated with these bubbles has a significantly harder spec- 
trum (dN/dE ~ E~ 2 ) than the inverse Compton emis- 
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Fig. 9. — Spectral energy distribution of three diffuse gamma- 
ray templates. The disk-correlated emission (green dashed) ap- 
proximately traces the inverse Compton and bremsstrahlung com- 
ponents. The spectrum of the uniform emission (dotted brown 
line) includes the isotropic part of the extragalactic background 
and cosmic-ray contamination. We have no t include d the dust- 
correlated emission or Loop I template as in ISu et al.l (2010), be- 
cause the expected spectrum is softer and we do not expect a signifi- 
cant contribution from these components at E > 100 GeV. Vertical 
bars show the marginalized 68% confidence range derived from the 
parameter covariance matrix for the template coefficients in each 
energy bin. 

sion from electrons in the Galactic disk, or the gamma- 
rays produced by decay of 7r° from proton-ISM collisions. 
We note that the morphology of the resolved gamma-ray 
cusp has a different shape from the bubbles, and the 
bubble structure has been largely cancelled out and not 
visible in the residual maps (as shown in e.g. Figure 
[5|). Thus we conclude that the new gamma-ray cusp has 
no obvious connection with the bubbles. We will also 
demonstrate this more explicitly in Section [U 

4. ENERGY SPECTRUM OF THE GAMMA-RAY CUSP 

In ^3] we argued for the existence of an excess with a 
concentrated cuspy shape in the GC at E = 120 — 140 
GeV. This motivates a more careful investigation of the 
energy spectrum and spatial distribution of the emission. 

4.1. Template regression 

As in our previous work (e.g., Su et al. 2010h . we con- 
sider linear combinations of spatial templates and com- 
pute the likelihood that the gamma-ray maps are de- 
scribed by a given linear combination. We fit a coefficient 
for each emission component using a multi-linear regres- 
sion of simple templates, one energy bin at a time. This 
technique provides an estimate of the energy spectrum 
of each component using as few physical assumptions as 
possible. 



By combining results from 16 logarithmically-spaced 
energy bands from 85 GeV to 200 GeV, we determine the 
spectral energy distribution for each component. In each 
fit, we model the "conventional" emission using three 
simple templates: a Galactic disk model, the Fermi bub- 
bles, and a uniform background. The Galactic disk tem- 
plate has the functional form (esc \b\) — 1 in lati tude and 
is a G aussian (pi = 80°) in longitude, as in ISu et al 



(2010). The disk model mostly accounts for gamma rays 
from the Galactic plane including those produced by 
point sources, ISM emission (tt° and Bremsstrahlung) 
and inverse Compton scattering. The uniform back- 
ground template absorbs the isotropic background due to 
extra-galactic emission and misclassified charged particle 
contamination, including heavy nuclei at high energies. 
We show maps of these templates in Figure El 

For each set of model parameters, we compute the Pois- 
son log likelihood, 



ln£ = y^kjlnfij - Hi- ln(fei!), 



(1) 



where Hi is the model counts map (i.e., linear combina- 
tion of templates times exposure map) at pixel z, and k is 
the map of observed counts. The last term is a function 
only of the observed maps. We compute parameter errors 
in the Gaussian approximation by inverting the matrix of 
second partial derivatives of — In C to obtain the covari- 
ance matrix, and taking the square root of the diagonals. 
The la Gaussian error corresponds to Aln£ = 1/2. 

Template-correlated spectra for the 3-template fit are 
shown in Figure The energy spectra show no signif- 
icant deviation from a power law for any of the three 
components. This fact, together with the distinct spa- 
tial morphology of the gamma-ray cusp suggests the need 
to include a cusp template (shown in Figure [8]) in the 
model. Inclusion of this template barely alters the de- 
rived spectrum of the first three components, but yields 
significant coefficients for the cusp from 110 to 140 GeV 
(Figure [10] and Table [Q). The fact that the bubble co- 
efficients show no such bump indicates that the bubble 
structure is unrelated to the ~ 130 GeV excess. We find 
that the surface brightness of the center of the cusp is 
nearly two orders of magnitude greater than that of the 
Fermi bubble structure, but only over a limited energy 
range. 

In Figure [10l we have done the fit with two different 
spatial masks: one excludes \b\ < 1° to avoid contami- 
nation from the Galactic disk close to the plane and the 
Galactic center (including the Galactic "ridge"); another 
one excludes much more of the Galactic disk with \b\ < 5° 
and |/| > 6° without masking the Galactic center region. 
We have obtained similar results for both data cuts. The 
cusp structure emits gamma rays with a luminosity of 
(3.2 ± 0.6) x 10 35 erg/s or (1.7 ± 0.4) x 10 36 photons/sec. 
The null hypothesis of zero intensity of the cusp compo- 
nent is ruled out by 5.0a. Given the energy resolution 
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Fig. 10. — Left panel: Spectral energy distributions of the templates listed in the figure legend. In the left panel, we use CLEAN events with 
|6| > 1° and all longitudes. Besides the disk-correlated emission (green), uniform emission (brown), and the Fermi bubble template (blue), 
the cusp component modeled as a FWHM = 4° Gaussian in the GC (red) has been included. Vertical bars show the marginalized 68% 
confidence range derived from the parameter covariance matrix for the template coefficients in each energy bin. Arrows indicate la upper 
limits. For reference, we overplot lines centered at 111 GeV and 129 GeV (dotted cyan) convolved with a three-Gaussian approximation of 
the LAT instrumental response (Edmonds 2011), and their sum (dotted black). The line centers and amplitudes are determined from a fit 
to the spectrum in the right panel (see text). Right panel: the same as the left panel but using data masking out \b\ < 5° and \l\ > 6°. 

unlikely statistical fluke, a dark matter annihilation sig- 
nal from the inner Galaxy is the most likely explanation. 
In another variant of the fit, we split the bubble template 
into two independent components in the fitting, high lat- 
itude (|6| > 30°) and low latitude (|6| < 30°). The pur- 
pose is to demonstrate that the low latitude bubble is 
also independent from the gamma-ray cusp. Again, we 
find no sign of a bump in the spectra of other diffuse 
gamma-ray components, but the cusp has a spectrum 
with an excess at 110 — 140 GeV and is consistent with 
zero in the other bins (Figure [TT]) . Instead of using CLEAN 
class, we have tried using SOURCE class for the likelihood 
analysis, and obtained similar results (Figure fl2|) . 

The energy spectrum of the cusp is consistent with 
a single spectral line (at energy 127.0 ± 2.0 GeV with 
x 2 = 4.48 for 4 d.o.f.). But a pair of lines at 110.8 ± 4.4 
GeV and 128.8 ±2.7 GeV provides a marginally better fit 
(with x 2 — 1-25 for 2 d.o.f.). We have compared the best 
fit one line and two line profile with the measured en- 
ergy spectrum in Figure [T3l The observation is compat- 
ible with a 140.8 ± 2.8 GeV WIMP annihilating through 
7Z and assuming rrih = 125 GeV (with x 2 = 3.33 
for 3 d.o.f.) or a 127.3 ± 2.7 GeV WIMP annihilating 
through 77 and 7Z ( with x 2 = 1-67 for 3 d.o.f.) (e.g., 



10- 7 



10- 8 L^ i_ 



j i i_ 



90 100 



120 140 160 
Photon Energy [GeV] 



-t — 



180 200 



Fig. 11. — Same as right panel of Figure [10] but splitting the 
bubble template into two regions one with |6| > 30° and the other 
with |6| < 30°. 



of the Fermi-LAT at E > 100 GeV, the spectral excess 
at 110 < E < 140 GeV is consistent with emission from 
one or two lines after cons idering the line-spread func- 
tion (LSF) (|Edmondsll20 1 1| ) . which strongly suggests the 
novel nature of the gamma-ray cusp as no known astro- 
physical process can produce this feature. Except for 
unexpected instrumental systematics or an increasingly 



Weiner Yavinl 120121 ). 



The gamma-ray cusp appears to possess a symmetric 
distribution around the Galactic center. To investigate 
whether there is any more extended cusp component con- 
tributing the excess at 120 — 140 GeV, we include an 
extra "outer ring" template as shown in Figure [HI The 
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Fig. 12. Same as Figure \W\ but using SOURCE events instead of CLEAN. The level of the uniform background is more than a factor of 2 
higher than that shown in Figure [TOl However, the resulting energy spectrum of the gamma-ray cusp (red dashed line) is quite similar to 
that shown in figures fTUI and ITTI 
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Fig. 13. — Same as right panel of Figure fTUI but overplot the 
best fit one gamma-ray line profile convolved with the instrument 
response, and compared with the best fit two-line profile. 

outer ring template is a FWHM=10° Gaussian with an 
8° radius hole in the center. Even with this freedom, 
there is no significant change in the cusp spectrum (Fig- 
ure fl4|) . There was no significant improvement of the 
likelihood for this model, and the spectrum of the outer 
ring is consistent with zero. Our conclusion is that the 
gamma-ray cusp is a distinct component, and is centrally 
concentrated. 

4.2. Trials factor 



Fig. 14. — Same as right panel of Figure [10] but including extra 
outer ring template. 

We use a trials factor of 300 for the single-line fits cen- 
tered on the Galactic center, and 6000 for fits that are off 
center. This choice is based on the fact that the LAT en- 
ergy resolution is ~ 10% over most of the energy range, 
and a line anywhere from 1 to 300 GeV would have been 
just as impressive, yielding 60 energy bins. Furthermore, 
a broader line (or two lines near each other) has an addi- 
tional trials factor. We allow an extra factor of 5, giving 
us 300. 

For fits that recenter the cusp in £, we would have 
considered any center with i\ < 5° interesting, and the 
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TABLE 1 

The template fitting coefficients and errors of the 
diffuse gamma-ray cusp correspond to the right panel of 
Figure [10] and right panel of Figure [12] The gamma-ray 
luminosity in each energy range is shown in the unit of 
keV cm - 2 s — 1 sr — 1 . 

centering must change by ±0.25° to make unit change in 
T5, so we take there to be 20 interesting spatial bins. For 
single-line, off-center fits we use a trials factor of 6000. 

The local significance of the centered Gaussian tem- 
plate is 5.0<j, obtained by summing the significance of 
each bin in quadrature (Table [1]). After diluting the p 
value corresponding to 5.0<j by a factor of 300, we obtain 
a significance corresponding to 3.7a. We note that, if the 
line is real, an additional 40% more data will be enough 
to obtain a 5.0<j detection, even with the trials factor of 
300. 

The previous trials factor is for a line anywhere with a 
range of widths. On the other hand, if one asserts that 
the line pair is from 77 and 7Z dark matter annihilation 
channels, the trials factor can be calculated as follows: 

The two 7 lines could have been anywhere between 
the Z mass and 300 GeV. There are 12 log-spaced bins 
of width 10% in that range. Also, there is no additional 
trials factor for a width; we simply have two lines, and 
we let their two amplitudes float as a free parameter, as 
well as the WIMP mass. There is also a factor from the 
number of line-producing scenarios: we could have seen 
a single line or two lines from 77 and 7Z or 7Z and 7I1. 
Considering these 3 scenarios, we assign a trials factor 
of 36 for the single-line fits. For off-center fits, we use 
36 x 20 = 720. 

5. DETAILED SPATIAL PROFILE OF THE CUSP 

The cusp template used in Section [4] is assumed to 
be centered on the GC, and the choice of 4° FWHM 
is somewhat arbitrary. We would like the data tell us 
what template to use, but there are so few photons, it is 
difficult to make sense of an unsmoothed map of counts. 
However, the smoothed maps indicate the cusp is slightly 
off center (to the W of the GC) and it is essential to follow 



this up. 

In this section, we consider individual photon events 
(not maps) and assume the exposure across the GC is 
slowly varying. We project the event locations into his- 
tograms of £ and b and study the distributions, finding 
parameters of a best-fit Gaussian. In this way, we may 
find the location, shape, and significance of the 130 GeV 
feature in a way that is independent of the previous sec- 
tions, though somewhat less principled because we do 
not explicitly use the exposure map. We will not use the 
results of this section to raise our claimed significance, 
but rather to emphasize that the cusp is centrally con- 
centrated, has a sharp spectrum, and is somewhat off 
center. 

We model the background spectrum to be dN/dE oc 
E~ 2 - 6 , with the amplitude in each bin set by the 10-50 
GeV average. The exposure map is a weak function of 
position and energy, and we neglect that variation in this 
analysis. An index of —2.5 gives a significantly worse fit 
by overestimating counts at 100 to 200 GeV. The 7r° 
emission is closer to —2.7 at lower energies, but —2.6 is a 
conservative choice, because assuming lower background 
at 130 GeV would make any excess more significant. 

We consider two spatial projections of the photon dis- 
tribution in the inner Galaxy. For the longitude projec- 
tion, we project the region \£\ < 10°, \b\ < 5°, in 0.5° 
bins of £ yielding a nearly flat distribution (blue line in 
top panels of Figure [T5|) . For the latitude distribution, 
we project the region — 5 < £ < 2, |6| < 10, in 0.5° bins 
of 6, finding the emission near the plane dominates (bot- 
tom panels of Figure [15]). From this we immediately see 
that the Galactic plane is much brighter than elsewhere, 
but the Galactic center is not particularly brighter than 
elsewhere in the plane at 10 to 50 GeV. 

In order to test for the existence of a bump, we com- 
pare the In Co for the null hypothesis to a model with an 
additional Gaussian of FWHM Fi, centered at £q with 
peak height Ag . We compute A In C = m(£/>Co) and ex - 
press results using the test statistic ( Mattox et al .11 199(1 ). 
TS = 2Aln£. The test statistic plays the role A% 2 
would play in a Gaussian problem. 

In comparing the ln£ of two models, one must account 
for the fact that the two models reside in different pa- 
rameter spaces. In our case, the null model space is a 
subspace of the other with 3 fewer parameters, obtained 
by setting Ai to zero. In this case, the TS distribution 
is simply the \ 2 distribution for 3 degrees of freedom. 

Using photons from all incidence angles, the addition 
of a Gaussian improves the TS by 36. This is not a 
6a result because of the 3 additional degrees of freedom. 
Rather, the probability that TS would be 36 or higher is 
p = 7.5 x 10 -8 , corresponding to 5.25a local significance 
(not including the global trials factor). The parameters 
of the Gaussian are Fg = 1. 4+0.4, A) = — 1-5 ± 0-3, and 
an amplitude corresponding to 14.0 photons (red line, 




Fig. 15. — Profiles for both £ and b. Even though the high-incidence-angle photons (0 > 40°; right) panels have half the exposure (9.7% 
vs. 19% for the left panels), they have more than half of the photons, and nearly the same TS due to lower off-line background leaking in. 
This demonstrates the statistical power of the high-incidence photons for line detection. See section [5] for a discussion of the significance. 



Figure [15]) . 

In the latitude direction, the fit is complicated by the 
concentration of conventional continuum emission in the 
plane. The cusp is not significantly offset in the b direc- 
tion, but sits in the region of highest background, so ad- 
dition of the cusp is not demanded as strongly by the fit. 
We introduce two new degrees of freedom, the amplitude 
and FWHM of a Gaussian centered at bo = 0. This yields 
TS = 28 .4 and p = 6.8 x 1(T 7 , corresponding to 4.8a (lo- 
cal significance). The maximum likelihood parameters of 
the Gaussian are = 3.9^o*7 an d A^ corresponding to 
16.1 photons. Both the t and b fits are roughly compati- 
ble with FWHM=3°, but there is a slight preference for 
an elongation of the cusp in the b direction. A careful 
study of this will require much more data. 

In Figure [15] (right panels) we also display the same 
plots for the high-incidence sample (0 > 40°). See Figure 
ITGl for such plots in 30 energy bins. The high-incidence- 
angle subsample contains half of the exposure time (9.7% 
vs. 19%) but due to better energy resolution (AE/E ~ 
0.06) has less background on the line, and therefore yields 
a TS almost as large as the full data. In this sense, most 



of the TS results from high events. This subsample 
would have yielded TS = 32.6(p = 3.9 x 10~ 7 , 4.93a) for 
the i profile, and TS = 26.1(p = 2.2 x lO" 6 , 4.59a) for the 
b profile. Although these are slightly worse p values than 
for the full data, they may actually be more persuasive 
due to the lower background. 

The fact that the cusp appears to be significantly off 
center implies that our spectral fit in the previous section 
erred by using a centered cusp template. In Figure [17] 
we show the measured energy spectrum of a 3° FWHM 
cusp template, centered at i = —1.5° and 6 = 0°. The 
local significance of this fit is 5.5a relative to the null 
hypothesis of zero intensity. This improvement is heart- 
ening; however, because of the extra parameter, the trials 
factor is now larger, diluting the significance. 

6. VALIDATION TESTS 

6.1. Assessment of line profile 

In section [3] we investigated the cusp emission by ana- 
lyzing maps in various energy bins. This allowed a sepa- 
ration of spectral components by morphology, but relied 
on an arbitrary choice of binning. The result - that there 
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Fig. 16. — Profile of high-incidence photon (6 > 40deg) longitude distribution for \b\ < 5°. The 0.5° bins have been smoothed by a 
3-bin box, and rescaled to arbitrary units of E 2 - 6 dN/dE, making the background disk emission constant with E. The background (blue) 
is the average in these units for 10 < E < 50 GeV. In each panel, the (non-negative) amplitude of a FWHM^ = 1.4° Gaussian centered at 
t = —1.5 is fit by maximizing the Poisson likelihood. The corresponding number of photons and test statistic (TS) are displayed. The only 
energy bin with significant emission is the 124.7-133.4 bin, centered on 129 GeV. See text for discussion of significance. The bin centered 
on 113 GeV is not significant by itself, but is compatible with a line strength of 1/3 to 1/2 that of the putative 129 GeV line. 



is a cusp of emission in the inner Galaxy - motivates an 
unbinned analysis of this region. 

In an unbinned analysis, one dispenses with arbitrary 
binning choices (size and shift) and instead analyzes in- 
dividual photon events. For example, the parameters of 
a well defined model may be estimated with no binning 
in space or energy. In the absence of a principled model, 
a compromise technique is to convolve a finely binned en- 



ergy histogram with some kernel and compare profiles of 
prospective lines with those expected for a true line, i.e. 
the instrumental response convolved with the smoothing 
kernel. 

In the case of LAT data this allows us to do an in- 
teresting reality check. Energy resolution of events at 
high incidence angle (6 ~ 60°) is a factor of ~ 2 better 
than that of normal- incidence photons, motivating the 
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Fig. 17. — Same as right panel of Figure [To] but with the cusp 
template centered at £ = — 1.5° and 6 = 0°. 
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Fig. 18. — Spectrum of emission within 4° of the cusp center 
(£,b) = (—1.5,0), excluding \b\ < 0.5°. High-incidence angle events 
(upper panel) have a factor of ~ 2 better energy resolution than 
those that enter the LAT close to normal incidence (middle panel) 
or the whole sample (lower panel). All three spectra have been 
smoothed by a Gaussian of 0.06 FWHM in AE/E, similar to the 
expected resolution of the upper panel. The continuum model is 
dN/dE ~ E~ 2Q , normalized at 20 < E < 50 GeV (blue dashed). 



following test. 

We select low incidence (6 < 30°) and high incidence 
(0 > 40°) photon samples. We restrict to those near the 
cusp center at (£,b) = (—1.5,0) (tjj c < 4°) but not in the 
plane (\b\ < 0.5°). We then convolve each with a ker- 
nel and compare them (Figure [T8|) . We adopt an LSF 
with a FWHM of A E/E 0.06 for high incidence and 
0.12 for low incidence ( Edmondsll2011 ), and in both cases 
convolve with another FWHM 0.06 Gaussian. After con- 
volution, the LSF is FWHM 0.085 for high incidence and 
0.134 for low incidence. Normalized Gaussians of these 
widths are shown for reference, normalized to the ex- 
pected line strength at 130 GeV. Maps constructed using 
only high incidence events are shown in Figure [191 

Note that: 

• The 129 GeV feature shape is strikingly similar to 
that expected for a line. The 111 GeV feature is 
unconvincing, but is also compatible with a line. 

• In some cases, fluctuations appear, but are not 
present in both low and high incidence spectra. 

This analysis did introduce some additional parameters, 
but we have made natural choices for them: The 68% 
containment radius of the cusp is approximately 4°, the 
Galactic ridge is about 0.5° thick, and the AE/E = 0.06 
smoothing kernel is similar to the LSF of the LAT at high 
incidence. Smoothing a spectrum by its LSF is often a 
good compromise between resolution and noise suppres- 
sion in the high- noise limit. Because these parameters 
are all fixed to natural values, there is no significant tri- 
als factor for this test, apart from the obvious one, that 
the lines could have appeared anywhere (Section [42]) . 

This test did not have to succeed. The fact that the 
high-incidence photon sample has sharper spectral fea- 
tures is important; if the high-# and low-# spectra in 
Figure [18] had been reversed, it would have been devas- 
tating for the line hypothesis. 

6.2. Null test: Galactic plane spectrum 

To emphasize that the line feature in Figure fT8l appears 
near the Galactic center and not elsewhere, we perform 
the same analysis on the Galactic plane (\b\ < 2°) away 
from the GC (^ > 5°). We find no indication of a line 
in either high-incidence or low-incidence photons (Figure 

Eq]). 

6.3. Null test: Earth emission photons 

Another null test is provided by the Earth emission 
photons. Cosmic-ray induced cascades in the Earth's 
atmosphere shower photons on the LAT at high zenith 
angle (Z > 108°). These provide another null test, as 
there is no reason for there to be a 130 GeV feature in 
the Earth emission spectrum. On average, no feature 
is seen (Figure [21]). However, there is a hint of a line 
at 130 GeV in the low-incidence events and one at 111 
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Fig. 19. — Same as Figure but using only events with high incidence angle > 40° which has better energy resolution. The maps have 
been smoothed to 5° (left two panels) and 10° (right two panels), and use CLEAN event class (upper two panels) and SOURCE event class 
(lower two panels). The resolved gamma-ray cusp structure is centered 1.5° W of the Galactic center. 



2000 



Galactic plane excl. GC (\|/ > 5°, |b| < 2.0° 




100 

E [GeV] 



2000 



GC spectrum (0 < 30°) 
111,129 GeV 
Expected LSF 




100 
E [GeV] 

Fig. 20. — Like Figure [18] except for |6| < 2° excluding photons 
within 5° of the Galactic center. There is no indication of a spectral 
feature near 130 GeV in the Galactic plane away from the Galactic 
center. 

GeV in the high-incidence events. These features also 
require a trials factor, but if they turn out to be robust, 
then serious concerns would be raised about instrumental 
artifacts giving rise to the observed lines. 

7. FITTING WITH DARK MATTER PROFILE 

In previous sections, we have modeled the line emission 
cusp with a Gaussian in an attempt to keep the analysis 
generic. However, WIMP annihilation is the most likely 
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Fig. 21. — Photons from cascades in the Earth's atmosphere 
(sometimes incorrectly called "albedo" photons), and with 10° of 
the GC, show a slight excess at 111 and 130 GeV also. Because 
these photons arrive at high zenith angle (Z > 110°), they tend 
to have a high incidence angle (median 6 = 63.2°). The low-# 
photons show a small bump at 130 GeV, and the high-# photons 
show a small bump at 111 GeV. The cuts were chosen to maximize 
these features, so interpretation of this plot requires a modest trials 
factor. 

explanation for gamma-ray line emission, and we would 
like to know if the signal expected from a reasonable 
DM profile fits the data better than a Gaussian. In this 
section, we repeat our template analysis with an Einasto 
profile as well as a set of modified NFW profiles. We also 
estimate the dark matter annihilation rate (<jv) required 
for each model. 

The gamma-ray intensity in a given direction is the 
line-of sight integral of the dark matter number density, 
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Fig. 22. — Same as right panel of Figure fTOl but using modified NFW profiles with inner power law index a for the DM density. The 
density is squared and projected along the line of sight to produce the cusp template for, from upper left to bottom right, a = 0.9, 1.0, 1.2, 
and 1.3, respectively. The fit prefers a = 1.2, yielding a 6.5a significance by combining the central 7 bins with > la each. 

earlier work, but a h igher value pQ ~ 0.4 GeV cm -3 
([Catena fc Ulliol [2QQ9h would reduce our (av) values by 
( 2 ) a factor of 1.7. 

For halo profiles p(r), we consider the Navarro- Frenk- 



squared, along a given direction, 



d<& 1 
dK 



with: 



J = 




(3) 



White (NFW) profile (jNavarro et al.H l 

Ps 



Pnfw(0 



(4) 



where the integral is over the cusp structure, {crv) 11 is 
the partial annihilation cross-section for XX ~~ ^ 77? m x 
is the WIMP mass (E 1 = m v ), ~ 8.5 kpc is the dis- 
tance from the Sun to the GC (|Ghez et al. 2QQ8| ). p(r) is 
the WIMP halo profile, p e ~ 0.3 GeV c m" 3 is the often- 
used WIMP density at the Solar system ([jungman et al.l 
19961 ). r = (s 2 + R% - 2s R Q cos £ cos 6) x / 2 is the Galac- 
tocentric distance and s is the line of sight distance. We 
use po ~ 0.3 GeV cm -3 to facilitate comparison with 



(r/r s ) a (l + r/r s ) 3 - a 
with r s = 20 kpc. The value of p s is determined by 
requiring p(r ) =0.3 GeV cm -3 . In Figure [22], we show 
the energy spectrum of the cusp structure using the NFW 
profile with various values of a. Interestingly, we obtain 
a 6.5a (5.6/5.9<j for one/two line case after trials factor 
correction) detection of the cusp assuming a = 1.2. This 
corresponds to a p value 10 4 times smaller than 5<r. 
We also consider the case of the Einasto profile, 

PEin as to(r) = p s exp{-(2/a £ ) [(r/r s ) aE - 1]} (5) 
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Models 


Before trials 


After trials (one line) 


Trials factor (one line) 


Gaussian (centered) 


5.0cr 


3.7a 


300 


Gaussian (off center, 9 > 40°) 


5.5(7 


3.7a 


6000 


unbinned £ 


5.2cr 


3.2a 


6000 


unbinned £ (0 > 40°) 


4.9a 


2.8a 


6000 


unbinned b 


4.8a 


3.5a 


300 


unbinned b (0 > 40°) 


4.6a 


3.2a 


300 


NFW a = 1.0 (off center) 


6.1a 


4.5a 


6000 


NFW a = 1.2 (off center) 


6.5a 


5.0a 


6000 


NFW a = 1.3 (off center) 


6.0a 


4.4a 


6000 


NFW a = 1.4 (off center) 


5.6a 


3.8a 


6000 


NFW a = 1.5 (off center) 


5.2a 


3.2a 


6000 


Einasto (off center) 


6.6a 


5.1a 


6000 



TABLE 2 

The detection significance of the gamma-ray cusp for various models. See £|4.2I for a discussion of trials factors. 



Models 


After trials (two line) 


Trials factor (two line) 


Gaussian (centered) 


4.3a 


36 


Gaussian (off center, > 40°) 


4.2a 


720 


NFW a = 1.0 (off center) 


4.9a 


720 


NFW a = 1.2 (off center) 


5.4a 


720 


NFW a = 1.3 (off center) 


4.8a 


720 


NFW a = 1.4 (off center) 


4.3a 


720 


NFW a = 1.5 (off center) 


3.8a 


720 


Einasto (off center) 


5.5a 


720 



TABLE 3 

The detection significance of the gamma-ray cusp structure with different models. See £14.21 for a discussion of trials 
factors. Masking out 0.5° around the GC area does not affect the results. 

8. HOW TO CONFIRM THIS SIGNAL QUICKLY: A 
MODIFIED SURVEY STRATEGY 

The fact that high-incidence-angle photons are supe- 
rior for line detection raises an exciting possibility: the 
scan strategy of Fermi-LAT could be altered for 1 year 
to confirm the 130 GeV line (if real) at 5<r with no trials 
factor, a significance that would be widely regarded as a 
discovery. If more than one line is present, the additional 
data would help characterize it. 

Fermi has usually scanned the sky in survey mode, 
observing the full sky every 2 orbits with occasional slews 
to targets of opportunity. This strategy is excellent for 
uniformity of full-sky coverage, but is far from optimal 
for collecting high-incidence- angle photons from the GC. 

From the spacecraft data files for the first 3.7 years 
(through week 202), it is a straightforward exercise to 
derive the fraction of the time the GC is accessible. We 
impose the constraints that the roll angle be within 35° 
of zenith (as in standard survey mode) and that the 
GC have an incidence angle of 45° < < 55°. We 
consider only times when the spacecraft is not in the 
South Atlantic Anomaly (SAA) and the DATA.QUAL flag 
is good. These constraints allow the GC to be observed 
40.6% of the time. This exceeds the exposure time of 
our (40° < < 60°) sample (observed 9.7% of the time) 
by more than a factor of 4. In other words, LAT could 
gather high photons from the GC 4x faster than it usu- 
ally does with a simple change to the observing strategy. 
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Fig. 23. — Same as right panel of Figure flOl but using an Einasto 
(a = 0.17) dark matter profile template to fit the data, we obtain a 
6.5a detection of the cusp only using the central 7 bins with > la 
each. 



with 



= 20 kpc and a = 0.17 (Einasto 1965 



Navarro et al.ll2004f ). This template provides a slightly 



better fit, with 6.6a (5.7/6.0<j for one/two line case af- 
ter trials factor correction) detection of the cusp. Table 
[2] and Table [3] lists the resulting significance for all the 
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Because we wish to point ~ 50° away from the GC any- 
way, it would still be possible to maintain coverage of the 
inner 110° radius about the GC. On the other half of the 
sky, survey mode would continue as usual. Although this 
change may be sub-optimal for some science projects, it 
would be a net benefit to the entire inner Galaxy, home 
to a great many scientific objectives. 

We believe the trade would be worth it. After 1 year 
of altered observing, we would have a sample of high 
incidence photons equal to the current sample, and could 
evaluate their significance directly, in the absence of any 
trials factor. In addition to confirming the line at 130 
GeV, other lines (113 GeV) may also become significant 
with additional data. At this point in the LAT mission, 
there are few ways for the instrument to dramatically 
improve its sensitivity to new physics in its remaining 
lifetime. This is a golden opportunity. 

9. DISCUSSION AND CONCLUSION 

Morphology: Using 3.7 years of Fermi-L AT data, we 
have identified a resolved gamma-ray cusp structure to- 
ward the Galactic center region. To reveal this struc- 
ture in a simple, model-independent way, we first take 
a linear combination of smoothed maps that cancels out 
continuum emission in the plane and inspect the resid- 
ual. We find that this structure only appears in the en- 
ergy range from ~ 120 GeV to ~ 140 GeV after search- 
ing 80 < E < 200 GeV maps (Figures EH and [T9l) . 
The FWHM of the cusp morphology, if modeled with a 
Gaussian, is < 4° and is u nrelated to the Fermi bu bble 
structure (as suggested by Profumo fc Lindenll2012 ). No 
other region of the sky reveals any significant excess in 
this energy range. 

Template regression: We perform a Poisson likelihood 
analysis to obtain the energy spectrum of the gamma- 
ray cusp. We create maps in energy bins 5.5% wide, 
and model each map as a linear combination of tem- 
plates including the cusp template, a uniform gamma-ray 
background, a thin gamma-ray disk tracing the Galactic 
plane, and the Fermi bubble structure. By modeling the 
data with a linear combination of templates and maxi- 
mizing the Poisson likelihood of observing the observed 
counts, we are able to separate the cusp emission from 
the Fermi bubbles and the Galactic disk. To the extent 
that the template is actually correct, this "matched fil- 
ter" gives an optimal estimate of the flux in each energy 
bin. The uncertainty estimates include marginalizing 
over the uncertainty in the other templates coefficients. 
We do not make a priori assumptions about the dark 
matter profile. For a template centered on the GC, we 
find that the cusp emits gamma rays with a luminosity of 
(3.2 ±0.6) x 10 35 erg/s, or (1.7 ±0.4) x 10 36 photons/sec. 
The null hypothesis of zero intensity is ruled out by 5.0a 
(3.7a with trials factor). 

Motivated by the apparent offset between the cusp and 
the GC, we repeat the template fit with a cusp template 



centered at (£,b) = (-1.5,0) (Figure [TTj), with the local 
significance rising to 5.5a, but with a larger trials factor, 
diluting the global significance. 

Line profile: The energy spectrum of the cusp struc- 
ture at E > 80 GeV is consistent with a single spectral 
line (at energy 127.0±2.0 GeV with x 2 = 4.48 for 4 d.o.f.) 



convo lved by the LAT energy response (jRajaraman et al 
2Q12h . A pair of lines at 110.8 ±4.4 GeV and 128.8 ±2.7 
GeV provides a marginally better fit (with x 2 — 1-25 for 
2 d.o.f.). These line energies suggest a WIMP annihi- 
lating to 77 and jZ. Fitting the two lines jointly, we 
find a WIMP mass of 127.3 ± 2.7 GeV with X 2 = 1-67 
for 3 degrees of freedom. The line pair is also compat- 
ible with a 141 GeV WIMP annihilating through 7Z 
and 7/1 f or rrih ~ 125 GeV, as in the "Higgs in Space" 
scenario ( Jackson et al. 2010h . We note that the uncer- 



tainty of overall absolution energy sca le of LAT is [-10, 
±5]% (|Fermi-LAT Collaboration! I2OI2I ). 



Interpretation as dark matter: Given the prop- 
erties of the gamma-ray cusp, and assuming it orig- 
inates from dark matter annihilation, constraints on 
the dark matter density profile and annihilation rate 
or cross section can be placed on various mod- 



els (e.g. Goodman et al. 2011: 


Jackson et al. 




2010: 


ICline 20121: Bertone et al. 2009: 


Dudas et al. 




2012; 


Gerinecer-Sameth & Koushiappasl 2012). At present, the 



extrapolation from the local dark matter density into the 
Galactic center is uncertain enough that there are large 
uncertainties in the cross section. We do not attempt to 
improve those estimates in this work. 

Comparison to previous work: The recent study by 



The Fermi-LAT collaboration: M. Ackermann et al. 
(|2012h put constraints on line emission from the inner 
Galaxy, excluding photons within \b\ < 5°. There- 
fore, their constraints are not in conflict with the 
signal claimed in this paper, which is mostly within 
5° of the GC. Ackermann et al. searched for dark 
matter gamma-ray line signal from 4.8 GeV to 264 
GeV (Fermi-LAT Collaboration: M. Ackermann et al 



2012h . updating the results fr om the earlier stud y with 



11 months Fermi-LAT data (|Abdo et all I2010D . The 
results from this study are obtained from two years of 
Fermi-LAT data and used the Pass 6 processing, which 
has somewhat worse background and instrumental 
systematics compared to Pass 7. 

The Galactic center is known to be a re- 
gion with diffuse high en ergy gamma-ray 



sion (|Aharonian et al.l l2006aT). especia lly the Galactic 



center ridge (|Aharonian et al.l l2006bf ). High energy 



gamma rays are produced by cosmic rays interacting 
with interstellar gas or giant molecular clouds in the 
central ~200 parsecs. We note that the resolved cusp is 
incompatible with a point source given that it extends 
to ~ 4°. To check for contribution from unresolved 
point sources, we have repeated our analyses masking 
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out regions with \b\ < 1°, and have obtained similar 
results, confirming that the cusp is not associated with 
the Galactic ridge. 

The future: The next version of Fermi-LAT data 
(Pass 8) will move us clo ser to realizing the full scien- 
tific potential of the LAT (|Atwoodll20ll ). The expected 
improvements include reduced backgrounds, increased ef- 
fective area, reduced point-spread function, better under- 
standing of the systematic uncertainties, and particularly 
extending the energy reach to higher photon energy. Fi- 
nally, the calorimeters of LAT have been demonstrated 
to effectively operate as a standalone instrument provid- 
ing imaging of the sky at E > 10 GeV, albeit with worse 
image resolution (at the level of 1°) with a large poten- 
tial in crease in the ef fective area at high energy and large 
angle ( Atwoodll2Q12| ). Such data would be valuable for 
future studies of the gamma ray line(s). 

A modified survey strategy: As discussed in section [HI 
the LAT could acquire significance on spectral line emis- 
sion from the Galactic center 4 times as fast as it cur- 
rently does with a change in observing strategy. In other 
words, a single year of observing the GC 40% of the time 
would yield more significance on these gamma-ray lines 
than all data currently in hand. A factor of 2 improve- 
ment would come from observing the GC more often, and 
another factor of 2 by maintaining a high incidence an- 
gle (45° < < 55°), yielding photon events with better 
energy reconstruction. We hope that the evidence pre- 
sented in this paper for a signal in general, and specifi- 
cally the superior value of the high-# events can be used 
to justify such a change. 

Data acquired from the new survey strategy, along with 
Pass 8 improvements, will be able to better constrain the 
morphology of the gamma-ray cusp and study in detail 
the line emission in 110 < E < 140 GeV, determine if 
there are two (or more) lines, and measure the ratios of 
the line strengths. Such measurements may be the key 
to differentiating between WIMP models and deepening 
our understanding of the dark sector. 

Is this a "discovery?": The local (without trials fac- 
tor) significance of our results is high: 5.0<j for the cen- 



tered Gaussian template fit to energy bins, 5.25cr for the 
unbinned Gaussian fit at {I = —1.5,6 = 0), and 6.5a 
by fitting the NFW dark matter profile centered at (- 
1.5, 0) with power index a = 1.2, with similar results 
for an Einasto profile. Given such high significances, it 
would be tempting to call this detection a "Discovery." 
However, we have a number of concerns: 

• The cusp is off-center by 1.5° (200 pc) and this was 
not expected. If the theoretical prior against this 
is strong, then this reduces confidence in the result 
- or at least a WIMP-related interpretation of the 
result. It is now an urgent question for simulators 
and theorists to determine whether this offset is 
unlikely or not in a spiral Galaxy with a significant 
bar. 

• The trials factor for this discovery is significant: 
The energy part is roughly 300 for a generic line 
search, though much less for the 77, 7Z scenario. 
The trials factor for the position offset could also 
be substantial. 

• Now that we have characterized what looks like a 
robust signal, the path to a truly convincing dis- 
covery is clear: modify the survey strategy to ac- 
cumulate signal on the GC as quickly as possible, 
and repeat the experiment with no trials factor. 
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